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LIBRARY ROCTINB K 10 » Igj 



Calculate Autocorrelation of a Tine Series 

Closed with four program parameters 

85 - total 

68 - main routine 

8 - constants at S3 

9 - Library Routine R 1 at SJ 

0,'1,2, -3 

-38 
Maximum error is 2 for the individual R. values. 

Approximately (3»5) (kN) milliseconds 

If a set of fractions or integers x,, x^, «.., x., •«.,. x 

are stored at locations y^ y 2 , ..., y^ •••7 n ** the form 

< \x ± \ < 10^ x 2* 59 "'- 



, this program computes 



"k- 



x i * X i+k " X i 



1+k 



<Ti .. G~i+k 



k -0, 1, 2, .,.,*' 
i = 1, 2, ..., N-k 



RESTRICTION 
REMARKS 



where ^ =J x ± 2 - xf 

and stores R^/p ^ locations z. > z p , •••, z , ..., z 

When the main routine is located at t the routine is 
entered by 

50 pF 



»• 



P+l 



26 tF 



with 



m' 



<KA 



y_ in location 
m* in location 1 
N in location JS3 
z, in A 



If (R1J/2 is not the form of storage desired, the values 
dRj^ may he stored by changing word 6l relative from 
S5 F JtOF to 75 68% ho F and storing the value 2d- as a 
fraction in relative location 68. 



DAT E ,7/13/35 RT: 9/13/98 
CODED B Y L. Aagenstein 
APPROVED BY J. P. Nash 



LOCAffKW 
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5 
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8 
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11 
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15 
16 



GREEK 

42 61L 
F4 IP 
42 6S3 
K5 F 

42 67L 
L5 F 
1*2 11 
42 8L 
42 10L 
L4 7S3 
42 5S3 
lH 3F 
41 1S3 
4l 2S3 
4l 3S3 
41 4S3 
41 S3 
L5 F 
40 F 
L4 S3 
40 S3 
L5 F 
40 IF 
L4 1S3 
40 1S3 
50 F 
75 IF 
S5 F 
L4 2S3 
40 2S3 
50 F 
75 F 
L5 3S3 
S4 F 



NOTES 
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K 10 



Set subroutine parameters 



Put x. in location and form Z x. 



Form Z x ± • x i+k 



Form Z x' 



LOCAfflCW 


ORDER 


iransB PAGE 2 


17 


to 3S3 
50 17 






18 


75 ar 

S5 F 




Ftallx^ 


19 


L4 4S3 
40 433 






20 


F5 8L 

42 8L 




Btep and test i and i+k 


21 


Y5 10L 
42 10L 






22 


LO 5S3 
32 8L 




i 


23 


L3 3F 
32 24L 




Test for k = 


24 


22 32L 
L5 733 


• 




25 


4o F 






26 


50 3F 
J5 3F 
1*2 3F 




Scale H to N x 2 211 > 4 2 x ± 2 


27 


L5 F 
00 2F 






28 


kO F 








LO 3S3 






29 


32 30L 
; 26 26t ' 




. 


30 


00 F 

F5 3F 




Set n shifts 


31 


k2 34L 
42 33L 






32 


41 3F 

L5 7S3 






33 


LO 3F 
00 F 




(H - k) (2 n ) In location 


3* 


40 F 




(N - k) (2 221 ) In location 1 




00 F 




r- 



K ; l# 



LOCATION 


ORDER 
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35 


40 IF 
50 S3 




(T^a^mss 




36 


85 F 
66 F 








37 


S5 F 
kO S3 








38 


50 1S3 
S5 F 








39 


66 F 
S5 F 




. (x 1+k ) 2" n in 1S3 


i 


ko 


ho 1S3 
50 2S3 
S5 F 






i 

i 
i 


hi 




(x i ' ^i+k^ 2 ^ 2S ? 


j 




66 IF 






! 
| 


I k2 


S5 F 
Uo 2S3 


• 




! 
i 
1 


M 


50 3S3 
S5 F 








hh 


66 IF 

S5 F 

^0 5S3 

50 hS3 




(7 X ) (2" 2n ) i* 3S5 
(^ +k ) (2" 2n ) in kS3 




h6 


S5 F 
66 IF 








hi 


35 F 

1*0 ^S3 








kS 


50 S3 
71 1S3 






' 


h9 


IA 233 
iiO 2S3 




Numerator x 2* iji 2S3 


j 


50 


50 S3 
71 S3 






j 


51 


I* 3S3 
00 IF 




<T ± in S3 


! 


52 


22 52L 
50 52L 




• 


{ 
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LOCATIQI 


ORDER 
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! K 10 


53 


26 SJ 
1*0 S3 






1 
| 
1 


5* 


50 1S3 
71 1S3 




• 




55 


iA kS3 
00 IF 




( ^l+k ) ( 2 ^ I/2 )^A 




56 


22 56L 
50 56L 








57 


26 SJ 
50 S3 








58 


4b IF 
75 IF 




2<T 2 (2" 211 ) in locatio* 1 




59 


kO IF 
50 2S3 








60 


S5 F 
66 IF 


' 


Store H / in 3k 




6l 


S5 F 
kO F 








62 


F5 3F 
k2 3F 




Step k 




63 


L5 IL 
\2 8L 




! 




6k 


lA- 3F 
42 1GL 




Reset loops and test 




65 


F5 6lL 
42 6lL 








66 


LO 6S3 
36 6L 








67 


22 67L 
22 F 
00 (S3)K 
00 F 
00 F 
00 F 
00 F 




Exit 
— 1 

1 





LOCATION 



CHHER 

00 F 

00 t 
00 J 
00 F 
00 F 
00 F 
HO S3 
L5 F 
35F 
J*0 F 
00 F 
00 F 

00 (SJ)K 
Routine R 1 - 116 
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HIGE 5 



K 10 



- Constant! and storage 
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APPROVED BT 




